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1. INTRODUCTION 

Spacebome and airborne radars are limited with 
respect to the mass and size of the instrument and the 
power available to operate it. As a consequence, dual- 
wavelength radars that require separate antennas and 
power amplifiers are expensive and often impractical. 
However, if the frequency difference can be reduced so 
that a single antenna and the same radio-frequency 
subsystem can be used for both frequencies, dual- 
wavelength Doppler measurements can be made with a 
radar of about the same size and mass as its single- 
frequency counterpart. 

In the first part of the paper we present calculations 
of the reflectivity factor differences as functions of the 
center frequency from 10 to 35 GHz and for frequency 
differences between -10% and 10% of the center 
frequency. The results indicate that differential- 
frequency operation at Ka-band frequencies (26.5 - 40 
GHz) provides relatively strong differential signals if the 
frequencies can be separated by at least 5%. Unlike 
lower frequency operation, the differential signals at Ka- 
band (both reflectivity and Doppler) are directly related 
to the median mass diameter. An important feature of 
the differential mean Doppler is that it depends only on 
the drop-size dependent part of the radial velocity. In 
principle, the mean and mean differential Doppler data 
from a nadir-looking platform can be used to infer 
vertical air motion and characteristics of the particle size 
distribution (Meneghini et al . , 2001). 

To test the instrument concept, the ER-2 Doppler 
radar (Heymsfield et al., 1996) was modified for 
differential frequency operation. Measurements by the 
modified radar, operating at frequencies of 9.1 GHz and 
10 GHz, were made using an 8° zenith-pointing offset 
parabolic antenna. Simultaneous data were taken with 
an optical rain gauge and an impact disdrometer. 
Measured and DSD-estimated values of the differential 
dBZ mean Doppler are presented in section 3. 


2. DIFFERENTIAL FREQUENCY CALCULATIONS 

For the results shown here, the hydrometeors are 
taken to be spherical and the scattering cross sections 
are calculated from Mie theory. We define the 
equivalent reflectivity factor difference, 5Z e (f, 8f), or 
more simply, the differential reflectivity, by 

8Za(f, 8f) = sgn(Sf) [dBZ e (f) - dBZ e (f+8f)] (1) 
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where sgn(x) = 1 for x>0 and -1 for x<0 so that 5Z e is 
always taken as the difference between dBZ e at the 
lower frequency to that at the higher frequency. The 
equivalent reflectivity factor at frequency f is given by 
dBZe(f) = 10 log io Z e (f) where 

Z e (f) =X 4 /[n S |Kw| 2 ] \ a b (A, D) N(D) dD (2) 

where a&(X, D) is the backscattering cross section 
(mm 2 ) of a spherical particle of diameter D at 
wavelength X and where N(D) is the drop size 
distribution (mm' 1 m' 3 ). For the calculations presented, 
N(D) is assumed to be an exponential. Because 5Z e is 
proportional to the log of a ratio of Z e measurements, it 
is independent of number concentration but dependent 
on the median mass diameter, D 0 . The differential 
mean Doppler, 8v(f, 5f), can be defined in a similar way: 

5v(f, 8f) = sgn(8f) [<v(f)> - <v(f+Sf)>] (3) 

where <v(f)> is the mean Doppler velocity at the radar 
frequency f. Shown in Fig. 1 are contour plots of 5Z e 
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Fig. 1 : Contour plots of SZ 9 in the f-6f space for rain with 
Do=2 mm (top) and snow with D 0s =4.31 mm (bottom). 
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Fig. 3 Maximum singular value plot for open and 
nominal closed loop system. 



(b) Performance path only 


Fig. 5 Maximum singular value plot for closed- 
loop system under different actuator failures 
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Fig. 4 Variation of H 0 0 norm with 5% errors in 
natural frequency, 500 cases 



Time, sec. 

Fig. 6 Time response of pointing under distur- 
bance for four conditions: Operational, Act-1 Fail- 
ure, Act-2 Failure, and Act-3 Failure 
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note that there are 3 mechanisms at work. In rain, at X- 
band frequencies, 8Z e is negative, i.e., dBZ« increases 
locally with f. Generally speaking, larger Do values are 
associated with larger values of |SZ 0 | although it must be 
noted that |5Z*| attains a maximum at D 0 =1.8 mm and 
begins to decrease thereafter. As the signals propagate 
into the rain 5Za becomes progressively less negative 
because the attenuation at the higher frequency is 
larger than that at the lower. Because 8Z e is positive in 
dry snow, a change in the sign of 5Z e occurs as the 
signals transit the melting layer. 
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Fig. 4. Height-time plots of <v(9.1 GHz)> (top), 5v(9.1, 
0.9) (center) and height-profiles of 5v (bottom). 


Somewhat similar characteristics can be seen in the 
differential mean Doppler results pictured in the center 
and lower panels of Fig. 4. in contrast to 5Z e , 5v is 
insensitive to differential attenuation. Moreover, 5v is 
approximately zero in the snow region. 

To validate the measurement concept, it is important 
to show the relation between the differential signals and 
the drop size distributions. To do this, we have used the 
disdrometer data and Mie theory to calculate the 
expected values of dBZ e , 5Z e , <v> and 8v. An example 
of comparisons between the disdrometer-derived dBZ e 
and 5Z e and the corresponding EDOP measurements is 
shown in Fig. 5. A 60-s shift of the EDOP data was 
introduced to account for the fact that the EDOP data 
were taken at a 400m height. Although comparisons 



Fig. 5: Comparisons of EDOP measurements (X) and 
disdrometer-derived (solid lines) values of dBZ(9.1 GHz) 
(top) and 5Z e (9.1 GHz, 0.9 GHz) over a 24 min period. 

between the EDOP data and DSD-derived estimates of 
SZa are noisy, the correlation is relatively good. This is 
encouraging in the sense that the (9.1 GHz, 10 GHz) 
combination is far from optimum, suggesting that sets of 
frequencies at Ka-band, separated by 5% to 10%, 
should yield stable estimates of 5Z e and 5v. 

4. SUMMARY 

Computations of the differential-frequency radar 
reflectivity factor and mean Doppler suggest that useful 
information on the drop size distribution and vertical air 
motion are feasible at Ka-band frequencies where the 
signal levels are relatively strong and directly related to 
the median drop diameter. As such, the measurement 
concept has applications to airborne and spaceborne 
sensing of rain and cloud. Experimental tests of the 
concept were conducted using the X-band EDOP radar 
where it was shown that measurements of 6Z e correlate 
fairly well with estimates derived from measured drop 
size distributions. 
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for rain (top) and snow (bottom) cases. For example, 
the values of 5Z e for the frequency pair (25 GHz, 26.25 
GHz) can be found at the (x, y) coordinates (25 GHz, 
+5%) yielding 8Z«= 0.45 dB in the case of rain and SZ*= 
0.6 dB for snow. For a fixed center frequency we can 
examine the behavior of 5Z e in the 8f-D 0 space by 
contour plots of the type shown in Fig. 2. Comparisons 
of the 14 and 35 GHz results reveal two advantages of 

f = 1 4 GHz 
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Fig. 2: Contour plots of 5Z e in 5f -Do space for fixed 
center frequencies of 14 GHz (top) and 35 GHz(bottom). 

operating at the higher frequency: larger signal strength 
and a 1-to-1 relationship between 8Z e and Do. At 14 
GHz, in contrast, more than one value of Do is 
consistent with the 8Ze measurement over certain 
ranges. Contours of n, the number of independent 
samples required for the mean differential signal level to 
be greater than the standard deviation of the 
measurement, are also shown on the figures. 

3. EDOP DIFFERENTIAL MEASUREMENTS 

Differential frequency measurements were made 
using the ER-2 Doppler radar (EDOP) in a ground- 
based, zenith-looking configuration. To modify the 
EDOP for differential measurements, the local oscillator 
(LO) at 9.6 GHz was disconnected and replaced with 
two synthesized HP-83640 sweep generators. These 
oscillators were used both to generate the transmit 
pulse and to mix the received signals down to the 


intermediate frequency (IF). The LO frequencies were 
chosen to produce transmit frequencies of 9.1 and 10 
GHz, a separation that was found to be the practical 
limit based on the performance characteristics of the 
EDOP traveling wave tube amplifier (TWTA). It should 
be noted that the generation, transmission, and 
reception of the two frequencies were done 
simultaneously (Bidwell et al., 2000). 

On March 27, 2000, an impact disdrometer and 
optical rain gauge were placed next to the zenith- 
directed antenna. Approximately 3 hours of data were 
collected. Because the modified radar was 

uncalibrated, the radar return powers were converted 
into radar reflectivities, by an additive constant, using 
the drop size distribution data near the beginning of the 
measurement period. Height versus time plots of 
dBZ e (9.1 GHz) and SZ e (f=9.1 GHz, 8f = 0.9 GHz) are 
shown in the top and center panels of Fig. 3. Displayed 
in the .bottom panel are selected vertical profiles of 5Z e . 
Corresponding plots of <v(9.1 GHz)> and 8v{9.1, 0.9) 
are shown in Fig. 4. The data show that the rain was 
primarily stratiform with a well-defined bright-band at a 
height of approximately 2.2 km. To understand the 5Z e 
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Fig. 3. Height-time plots of dBZ e (9.1 GHz) (top), 5Ze(9.1 
GHz, 0.9 GHz) (center) and height-profiles of 8Z e 
(bottom). 


